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^^ \ Abstract: We propose a simple method to simulate W/Z + jet productions at hadron 

■^ ' collisions. The simulation can be done by using existing tools with some modifications 

o ■ allowed to users. W/Z + 1 jet events are generated using an ME-based event generator at 

X.' the tree level. The divergence at low px is suppressed by using the Sudakov form factor. 

D . PS is added with an appropriate consideration for PS-ME matching. The simulation for 

the VF + 1 jet production shows a smooth matching with the W^ + jet simulation at 

low pTS. The Z -\- 1 jet simulation in the Tevatron Run I condition well reproduces the 

rS ' experimental measurement on the px spectrum of Z. 
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1. Introduction 

The simulation of multiple hadron jet (multi-jet) production is one of the most serious 
problems in physics analyses at high-energy hadron collider experiments. Our experimen- 
tal reach extends to heavier objects as the collision energy increases. Once such heavy 
objects decay hadronically, they frequently produce multi-jet final states. Similar multi-jet 
configuration originating from well-separated emission of light quarks and/or gluons can be 
easily produced by non-resonant QCD interactions. They may blind or distort interesting 
heavy object signals. We have already encountered this problem in the study of hadronic 
top-quark decays at Tevatron. The problem will become much more severe for heavier 
objects (Higgs boson(s), SUSY particles, etc.) expected at forthcoming LHC experiments, 
since in principle the signal frequency decreases as the cm energy of the hard interaction 
increases while the QCD activity remains nearly constant. Understanding of the QCD 
multi-jet production will become crucial for the study of such heavy objects. Despite that, 
the technique is not well established for the simulation. 

The simulation of hadron collision interactions consists of two components: parton 
showers (PS) describing relatively soft regions in the initial-state and the final-state in- 
teractions, and a simulation of hard interactions based on perturbative matrix-element 
(ME) calculations. The former is a 3-dimensional model based on the factorization theory, 
adding the contributions of large collinear components to all orders. The matching be- 
tween the two components is important, but is not trivial since theoretical discussions are 
usually made only at the collinear limit where transverse behaviors are ignored. A certain 
model-based discussion is necessary to construct a consistent multi-jet simulation. 

Recently, the CKKW technique originally proposed for the final-state jet production 
at e+e~ collisions |||] has been extended to the initial state B, and implemented in an 



event generator for hadron collisions y]. This technique consistently adds tree-level ME 
calculations for a certain hard interaction process associated with jet, 1 jet, 2 jets, ..., to 
provide an exclusive event sample above a certain resolution scale. A PS is applied to a 
limited phase space which ME calculations do not cover, with a careful consideration for 
the matching. The concept is also implemented in different ways in other event generators 

i- 

CKKW provides us with an exclusive multi-jet sample. Despite, in many cases, users 

are interested only in an inclusive behavior of a certain number of leading jets. They ignore 

non-leading ones in order to make the evaluation as free from theoretical and experimental 

ambiguities in soft regions as possible. We expect that there must be a simpler method 

for such applications. In this paper we propose a method to obtain a simulation sample of 

inclusive n-jet events by using an n-jet ME event generator with the help of an appropriate 

PS simulation. 

In the present study we focus on the simplest case, the VF + 1 jet production. Namely, 

we are interested only in the leading-jet behavior produced in association with the W boson 

production. This process deals with the jet emission from the initial-state partons only. 

The discussions would be able to extend to the final state with appropriate replacements 

of the parameters and formulae. We use only those tools which are publicly available. 

Modifications are applied at the level where ordinary users are allowed. We apply our 

method to the Z + 1 jet production in the Tevatron Run I condition for a comparison 

with experimental measurements. The Z production is used only because the momentum 

measurement is expected to be less ambiguous than the W production. 

2. W + 1 jet simulation 

The VF + 1 jet production is simulated at the tree level by using GR@PPA version 2.76 ||5| 
together with the PS in PYTHIA version 6.212 |g in the present study. The W bosons 
always decay to the pair of an electron and a neutrino. GR@PPA includes this decay in 
the ME calculation. The sample program in the GR@PPA distribution package is used 
for interfacing GR@PPA and PYTHIA. Studies are done for the LHC condition, proton- 
proton collisions at the cm energy of 14 TeV. A CTEQ6L1 Q routine included in the 
sample program is used for PDF. 

Both the initial and final state PSs are activated in PYTHIA with the default set- 
ting, while the hadronization and the multiple interaction are deactivated for simplicity 
(MSTP(81) = and MSTP(lll) = O). The QED radiation is also deactivated (MSTJ(41) 
= l). A jet clustering (PYCELL in PYTHIA) is applied to the parton-level events, where 
the detector is assumed to cover the full azimuth ((/>) and the pseudorapidity (77) up to 4.5 
in the absolute value with a granularity of about 0.1 in both cj) and r]. The jets are recon- 
structed using a cone algorithm with the half-cone size of 0.4 in R and with the transverse 
energy (Et) threshold of 10 GeV. 

A naive simulation of this process has an obvious difficulty. In hadron collision simu- 
lations we have to give four energy scales: renormalization scale (/i/j), factorization scale 
(fj-p), and energy scales for the initial-state PS (/^isr) and the final-state PS (/ipsR)- It is 
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Figure 1: Naive simulation oiW + 1 jet production at LHC. The pt distribution of the W boson 
and the Et distribution of the leading jet are plotted. The M^ + 1 jet events are generated by 
GR@PPA 2.76 with the pr cut of 5 GeV/c, and the PYTHIA PS is apphed to the initial and final 
states with the default setting. A Et cut of 10 GeV is imposed in the jet reconstruction. Results 
for two extreme cases of the energy-scale choice are shown: ^ = px with solid circles, and /i — mw 
with open circles. Histograms show the results from the VF + jet simulation with fi = mw- 



usually said that they should be set to the typical energy scale of the interaction. However, 
there are two energy scales in the TV + 1 jet production: the VF-boson mass (mw) and the 
transverse momentum (pt) of the jet. They can be quite different from each other if we 
allow very low px (~ 10 GeV/c) for the jet. 

Figure |l| shows the simulation results for the px distribution of W and the Et distribu- 
tion of the leading jet (the highest Et jet). The ME events are generated with a minimum 
Pt cut of 5 GeV/c in the cm frame of the hard interaction. Distributions are shown for two 
extreme cases of the energy-scale choice. All energy scales are defined to be identical and 
set equal to the jet pT in the cm frame (pt) in one case, and the W mass (the invariant mass 
of the electron and neutrino) in the other. These choices are non-standard in GR@PPA. 
The user-define option (ICOUP = IFACT = 6) is selected, and the energy-scale parameters 
(GRCQ and GRCFAQ) are appropriately set in the subroutine GRCUSRSETQ. We frequently 
see similar simulations in previous W + jets analyses. The figure shows that the difference 
is significant not only in the absolute value but also in the shape. 

The prediction from the VF + jet event generator in GR@PPA 2.76 is also shown in 
Fig. ffl. The same PYTHIA PS is applied with all the energy scales set equal to the W 
mass. The finite pT of W and all jet activities are generated by PS in this case. The low 
Pt behavior must be well described by this simulation. However, neither the two VF + 1 
jet simulations agrees with it at low pts. The fi = pT simulation looks better around pT = 
20 GeV/c, but deviates at further low pts. This is unavoidable since the W + 1 jet cross 
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Figure 2: One of the Feynman diagrams for which the VF + 1 jet ME is evaluated. 

section diverges at pr = 0. 

Figure ^ illustrates one of the Feynman diagrams for which the T^ + 1 jet ME is 
evaluated. In the present inclusive analysis, the jet in the ME should be identical to the 
leading jet that we observe. This means that there should not be any jet having px larger 
than that. The naive simulation with /x = my/ is apparently inadequate from this point of 
view. The initial-state PS frequently generates high pT jets, resulting in large cross sections 
in medium pT and Ex regions in Fig. [l|. 

The excessive PS activities can be eliminated by the choice of /i = pt- This leads to a 
better behavior in low pT and Et regions. The divergence at further low p^s appears only 
because we stop the perturbation at a limited order, the lowest order in this case. Higher 
order corrections would produce higher pT jets and suppress the contribution of the TV + 1 
jet ME at low p^s. This is considered to be the mechanism to make the actual observation 
finite. The suppression can be evaluated to all orders of the coupling constant in a collinear 
approximation in terms of the Sudakov form factor. This is the way usually adopted in 
QCD calculations having multiple energy scales, and the way used in PS. Therefore, the 
application of such a suppression with an appropriately evaluated Sudakov form factor 
should modify the VF + 1 jet simulation to match the 1^ + jet simulation at low pts, and 
to be finite even at pT = 0. We look for a suitable expression of the Sudakov form factor 
in the next section. 

Although the choice oi ^ = pT looks better, it is too naive to construct a realistic 
simulation. We reconsider the choice of energy scales and discuss about the matching 
between PS and ME in later sections. 



3. Sudakov form factor 



The Sudakov form factor for quarks can be described as 



Sq{Qi,Q2) =exp 
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target parton 
Pj=(Ej,Pt,0,Pi^) p^ = (£2,0,0,- Pi^) 

z = s'ls __ 

^S' = ipi + Plf 




evolving parton \ radiation 

Pl ={E^, 0,0, Pi} A P4 = (E4-PT^0,Pi^-p'i^) 

S^iPl+Plf 
Figure 3: The definition of the PS sphtting kinematics that we adopt. 

where the splitting function for quarks with the radiation of a timelike gluon is 

P,^,{z) = Cf^^ (3.2) 



at the leading order, with Cp = 4/3. Equation ( |3.1| ) represents the probability that a quark 
or an anti-quark survives without any radiation from an evolution parameter value of Qi 
up to QI- However, the quantity is unphysical since the result depends on an artificial 
cutoff e. A smaller e gives a smaller Sq. This corresponds to a natural explanation that 
the non-radiation probability becomes smaller if we allow softer radiation, but the softness 
is not well defined. It is necessary to introduce a definition of the splitting kinematics in 
order to derive a physically meaningful quantity. 

Here we examine the definition in PYTHIA for the initial-state PS Q, where Q^ is 
exactly the virtual of evolving partons. Details of the definition are shown in Fig. ^. The 
four-momenta are so defined that pf, p\ and p\ < 0, and pi > 0. The splitting parameter 
z is the ratio of the squared invariant mass of the collision system after and before the 
splitting, given as 

S [Pi + P2) 

This definition preserves the relation Shard = X1X2S with x given by the product of all 
z values in each beam. The squared four-momenta, s and z are the inputs. The other 
parameters can be derived from them. 

In this definition, the transverse momentum of the splitting can be described as [0] 



p'T = Ei-p'^-pi (3.4) 



where 
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Figure 4: The integration area defined by the 6 function in Eq. ( |3.7| 



In order to simplify the relation we assume the radiation to be massless (p^ = 0), and 
take the limit p^—^Q since fundamental properties should not depend on not-well-defined 
target virtuality. Furthermore, we take the limit p^ 



pI 



-Q since Eq. ( |3.1| ) indicates 
that we can take the virtuality step as small as we want. These approximations lead us to 
a relation, 

(1-Z)2Q2 



Pt 



1 + ^ 



If we assume Q^ <^ s, this can be further simplified to 

PT = {1- z)Q. 



(3.5) 



(3.6) 



This is different from the usually quoted relation p^ = (1 — z)Q'^ , but similar to the one 
obtained for the HERWIG PS [|. 

If we assume the relation ( |3.6| ) , the Sudakov form factor giving the no-radiation prob- 
ability above a certain px cut (pTmin) can be calculated as 



SqipTmini Qr 
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"Q2" 1^2^ Pg^q{z)e{pT - PTn 



(3.7) 



The integration area defined by the 9 function is schematically illustrated in Fig. ^. The 
9 function naturally gives an upper bound of the z integration and an lower bound of the 
Q^ integration. Therefore, there is no need to introduce any artificial cutoff. Namely, the 
result is well-defined and must have a physical meaning. The upper bound Qmax should be 
given by the hardest energy scale in the considered interaction. It would be natural to take 
the VF-boson mass for non-hard radiation in the M^-l- 1 jet production. Although Eq. ( |3.7] ) 
is numerically calculable, it is rather complicated for the application to actual simulations. 
Let's start from the radiation function which the Sudakov form factor is based on. The 
radiation probability is usually expressed as 






(3.8) 



If we assume the relation (p.q), we obtain another expression that 



Using this expression, the Sudakov form factor having the same meaning as Eq. 
be written as 



SqipTmin, Qmax) = exp 
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(3.10) 



The upper bound of the z integration is given by the relation px = {I — z)Q < (1 — ^)<3max- 
Since it is natural to define a^ to be a function of px of the splitting, the z integration 



is easy to perform. Assuming the leading order function (3^), we obtain 

3-4e + e2 



Sq{pTmm,Qm£ix) = exp 



vr PT \ e 



PT- 



dpT 



(3.11) 



with e = pr/Qmsix- This leads to the expression used in the CKKW method [Q] if we take 
the limit e — > in the non-divergent term, although the definition of the parameters is 



slightly different. The difference between Eq. ( |3.7] ) and Eq. ( 3.11 ) is only in the definition 
of the integration variables. They give exactly identical answers. 



4. Suppressed W^ + 1 jet simulation 

For the test of the suppression, VF + 1 jet events are generated using GR@PPA 2.76 with 
the energy scale choice of /i/j = ^p = px, where px is the transverse momentum of the 
jet and W in the cm frame of the hard interaction. The choice /x/j = pT merely means 
that the coupling at the parton splitting is evaluated at this energy scale. A CKKW-like 
correction has to be applied when we have multiple QCD vertices. The other conditions 
are the same as those in the simulations leading to Fig. ffl. 
The weight factor for the suppression is defined as 



w 



Sq{pT,Qma.xy 



(4.1) 



with Eq. (3.11), where we use the first-order expression for the strong coupling: 

An 



asipx) 



with /3o = 11 — 2nf/3 (nj = 5) and A 
hard interaction scale is given by 



PoHPt/^^) 

0.0883 GeV. This gives as{mz) 



Q. 



mriW) = ^Jml^, + 



Pt- 



(4.2) 
0.118. The 

(4.3) 



This naturally connects reasonable definitions at two extreme cases: Qmax = ^^w ^^ small 
Pts and Qmax = Pt at large pts. The Sudakov form factor is squared since the no-radiation 
condition must be required to two incoming quarks converted to the W boson. Events 



generated according to the ordinary W + 1 jet cross section are accepted in proportion to 
this weight in the subroutine GRCUSRCUT of GR@PPA. Though this degrades the event 
generation efficiency, the degradation is moderate in GR@PPA because the same routine 
is used in the initiaUzation stage for optimizing the event generation. 

We have to apply PS also to these VF + 1 jet events in order to obtain realistic events. 
Here, an appropriate care is necessary to accomplish a reasonable matching between the 
PS and ME. In principle the PYTHIA PS is equivalent to the numerical evaluation of the 
integral in Eq. ( |3.lD . The Sudakov form factor is evaluated by carrying out the integration 
over the shaded area in Fig. §. Thus, this area is covered by ME if pTmin = Pt', i-e., no 
radiation there. PS has to cover the rest of the area. We can easily see that the matching 
cannot be achieved by a simple choice of the PS energy scale. A certain rejection method 
is necessary to apply, as is done in CKKW, when we use an ordinary Q^-ordered PS. 

We adopt the following method in the present study: the PYTHIA PS is applied with 
the energy-scale choice of fips = rhT{W). This can be done by explicitly setting the variable 
SCALUP of the LHA common HEPEUP to this value in the subroutine UPEVNT. After the PS 
is added we investigate the parton information in the PYJETS common of PYTHIA. The PS 
is re-applied to the same event if there is any radiation from the initial-state partons with 
Pt ^ Pt- This is an approximation since the examined px does not directly correspond to 
the Pt of each splitting. However, this method can be applied without any modification 
to PYTHIA routines and, what is more, it provides us with a reasonably good results as 
will be shown in the following. In any case, it is impossible to achieve a perfect matching 
using Q^-ordered PSs. There is no splitting to produce spacelike gluons in the VF + 1 jet 
ME, while in PS such a splitting may produce the highest pj- jet although the probability 
must be very small. This rejection does not significantly affect the event generation speed, 
because the PS simulation is fast and the average number of trials is only 1.3 in the present 
condition. 

Figure ^ shows the pt{W) and £^T(leading jet) distributions of the suppressed W + 
1 jet simulation. Only low to medium pj- and Et regions are shown in the figure. The 
suppression is not significant in higher px (Et) regions. The results from a naive simulation 
{fi = Pt) is overwritten with dashed histograms to show how the suppression works. We 
can see a smooth matching with the TV + jet simulation (histograms) at low pT (Et) 
around 15 GeV. This shows that our suppression method is reasonable. 

Two results are plotted in the figure. The filled circles show the results of the sim- 
ulation implementing the above PS-ME matching method, while the open circles show 
those without matching. The PS energy scale is set to px in the latter in order to avoid 
the overlap with ME. The threshold behavior is more strongly smeared in the matched 
simulation. 

The difference between the matched and unmatched simulations is more clearly seen 
in Fig. P, where pt{W) is plotted for jet-tagged events. The distribution below the Et 
threshold (10 GeV) shows the collective effect of the applied PS. The result for the VF + 
jet simulation obtained from the same analysis is overwritten in the figure. The matched 
simulation is in good agreement with the VF + jet simulation below the threshold. This 
means that our matching method is well performed; namely, the PS added to the VF + 1 
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Figure 5: Suppressed M^ + 1 jet simulation in the LHC condition. A pT cut of 5 GeV/c is applied 
in the hard-interaction generation, and the jet reconstruction is applied with an Ex threshold of 
10 GcV. Filled circles show the results of the simulation implementing the PS-ME matching, while 
open circles show those without matching but with fips = pt- Dashed histograms are the naive 
simulation results for the energy-scale choice of /x = pr- Solid histograms show the results from the 
VF -|- jet simulation with fi = mn/. 

jet simulation is quite similar to the non-leading PS in the VF -|- jet simulation. Since 
the leading behavior is already well matched, the suppressed W^ -|- 1 jet simulation is now 
indistinguishable from the W + jet simulation at low pxs. 

Figure 1^ shows the total cross section of jet-tagged events for a practical jet Et thresh- 
old (20 GeV), as a function of the px cut in GR@PPA. Observable quantities should not 
depend on artificial conditions in the simulation such as the pT cut. The result shows 
that the jet-tagged cross section is stable against the variation of the pT cut if it is set 
reasonably small. The pT cut of 10 GeV/c looks sufficient in the present study, though 
further smearing due to the hadronization and detector effects may alter it. 



5. Z production at Tevatron Run I 

In this section we apply our suppressed simulation to the Z production at Tevatron for a 
comparison with experimental data. The simulation for the Z production (66 < m{Z — > 
e'^e~) < 116 GeV/c^) is compared with CDF data at Tevatron Run I p!o| ] in Fig. |8|. 
The result from the suppressed Z + 1 jet simulation for the Tevatron Run I condition 
{pp collisions at 1.8 TeV in the cm energy) is shown with the histogram. A pT cut of 5 
GeV/c is applied again. The simulation is essentially the same as that for the VF -|- 1 jet 
production, except for the overall normalization. The simulation result is multiplied by a 
factor of 1.48, the ratio between the measured inclusive Z production cross section and the 
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p/W) (GeV/c) 

Figure 6: The pr distribution of W for 
jct-taggcd events. Filled circles show the re- 
sult of the simulation to which the PS-ME 
matching method is applied, while open cir- 
cles without the matching but with the PS 
energy-scale choice of p.ps ~ pr- The his- 
togram shows the result of the same analysis 
applied to the VF -I- jet simulation. 



E.p(leadingjet)>20GeV 



p^^,„ (GeV/c) 

Figure 7: Total cross section of jct-taggcd 
events {Et > 20 GeV) as a function of 
the pT-cnt in the hard-interaction genera- 
tion. Filled circles show the results from the 
suppressed W + 1 jet simulation with the 
PS-ME matching. Open circles show those 
we obtain if we omit the rejection procedure 
in the matching. 



data points: CDF, PRL 84, 845 (2000) 
histogram: simulation witii K-factor = 1.48 
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Figure 8: Comparison with CDF data at Tevatron Run I for the pT spectrum of Z. The histogram 
shows the simulation result for the Z + I jet production with the pT cut of 5 GeV/c. 



corresponding tree-level prediction by the Z + jet generator in GR@PPA. We can see a 
very good agreement from low (~ 15 GeV/c) to very high (~ 100 GeV/c) px regions. It 
should be noted that we have obtained this simulation result without any tuning. 

Unfortunately jet spectrum data are not available for the comparison. There would be 
a difficulty in the presentation of experimental data since jets are experiment and analysis 
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dependent. By the way, there is no reason that experimental data much differ from our 
simulation, since the agreement is very good in the pt{Z) spectrum. At the 0th order, 
the leading jet should be balanced with Z in px- Though higher order effects violate 
the balance, the dominant leading-order corrections are taken into account by applying 
a PS in our simulation. Thus, significant deviation can emerge only if we examine a 
quantity relevant to the next-to-leading jet. Note that the events generated by GR@PPA 
are passed to PYTHIA with the energy scale choice of fips = rhT{Z). The final-state PS 
is implemented up to this scale since no rejection is applied to it. 

6. Conclusion 

We have introduced a simple method for simulating inclusive W/Z + jet productions at 
hadron collisions. The divergence of the cross section at low px is suppressed by using the 
Sudakov form factor. PS is added to the events generated by an ME-based event generator 
at the tree level with an appropriate consideration for matching. The simulation can be 
done using existing tools, GR@PPA and PYTHIA. Necessary modifications can be applied 
at the level where ordinary users are allowed. 

The results from our VF -|- 1 jet simulation show a smooth matching at low pts with the 
Ty -|- jet simulation, in which jet activities are totally generated by PS. The simulation 
of Z + 1 jet production in the Tevatron Run I condition well reproduces the CDF data for 
Pt{Z) measurement. 

The present study deals with the jet (parton) radiation from initial-state partons only. 
An extension to the final state is necessary to apply our method to multi-jet productions. 

The method applied in the present study is similar to the inclusive treatment for the 
highest jet-multiplicity events in CKKW Q. The success of our method suggests that 
the inclusive simulation of n-jet events in CKKW might be enough to obtain a reasonable 
inclusive n-jet sample, if the resolution scale can be lowered to the level of the px cut in the 
present study. The contribution from other multiplicity events might become insignificant. 
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